The positivity rates and drug resistance patterns of Mycobacterium tuberculosis using nucleotide MALDI-TOF MS assay among suspected tuberculosis patients in Shandong, China: a multi-center prospective study.

Objective: To investigate the positivity rates and drug resistance characteristics of Mycobacterium tuberculosis (MTB) among suspected tuberculosis (TB) patients in Shandong Province, the second-largest population province in China. Methods: A prospective, multi-center study was conducted from April 2022 to June 2023. Pathogen and drug resistance were identified using nucleotide matrix-assisted laser desorption ionization time-of-flight mass spectrometry (nucleotide MALDI-TOF MS). Results: Of 940 suspected TB patients included in this study, 552 cases were found to be infected with MTB giving an overall positivity rate of 58.72%. Total of 346 cases were resistant to arbitrary anti-TB drug (62.68%), with Zibo (76.47%), Liaocheng and Weihai (both 69.23%) ranking top three and TB treatment history might be a related factor. Monoresistance was the most common pattern (33.53%), with isoniazid the highest at 12.43%, followed by rifampicin at 9.54%. Further analysis of gene mutations conferring resistance revealed diverse types with high heteroresistance rate found in multiple anti-TB drugs. Conclusion: A relatively high rate of MTB positivity and drug resistance was found in Shandong Province during and after the COVID-19 pandemic, indicating the need for strengthening rapid identification of species and drug resistance among suspected TB patients to guide better medication and minimize the occurrence of drug resistance.

Transmission dynamics and phylogeography of Mycobacterium tuberculosis in China based on whole-genome phylogenetic analysis.

OBJECTIVES: This study aimed to describe the lineage-specific transmissibility and epidemiological migration of Mycobacterium tuberculosis in China. METHODS: We curated a large set of whole-genome sequences from 3204 M. tuberculosis isolates, including thousands of newly sequenced genomes, and applied a series of metrics to compare the transmissibility of M. tuberculosis strains between lineages and sublineages. The countrywide transmission patterns of major lineages were explored. RESULTS: We found that lineage 2 (L2) was the most prevalent lineage in China (85.7%), with the major sublineage 2.2.1 (80.9%), followed by lineage 4 (L4) (13.8%), which comprises major sublineages 4.2 (1.5%), 4.4 (6.2%) and 4.5 (5.8%). We showed evidence for frequent cross-regional spread and large cluster formation of L2.2.1 strains, whereas L4 strains were relatively geographically restricted in China. Next, we applied a series of genomic indices to evaluate M. tuberculosis strain transmissibility and uncovered higher transmissibility of L2.2.1 compared with the L2.2.2 and L4 sublineages. Phylogeographic analysis showed that southern, eastern, and northern China were highly connected regions for countrywide L2.2.1 strain spread. CONCLUSIONS: The present study provides insights into the different transmission and migration patterns of the major M. tuberculosis lineages in China and highlights that transmissible L2.2.1 is a threat to tuberculosis control.

Epidemiology and Association Rules Analysis for Pulmonary Tuberculosis Cases with Extrapulmonary Tuberculosis from Age and Gender Perspective: A Large-Scale Retrospective Multicenter Observational Study in China.

Background: Tuberculosis (TB), a multisystemic disease with protean presentation, remains a major global health problem. Although concurrent pulmonary tuberculosis (PTB) and extrapulmonary tuberculosis (EPTB) cases are commonly observed clinically, knowledge regarding concurrent PTB-EPTB is limited. Here, a large-scale multicenter observational study conducted in China aimed to study the epidemiology of concurrent PTB-EPTB cases by diagnostically defining TB types and then implementing association rules analysis. Methods: The retrospective study was conducted at 21 hospitals in 15 provinces in China and included all inpatients with confirmed TB diagnoses admitted from Jan 2011 to Dec 2017. Association rules analysis was conducted for cases with concurrent PTB and various types of EPTB using the Apriori algorithm. Results: Evaluation of 438,979TB inpatients indicated PTB was the most commonly diagnosed (82.05%) followed by tuberculous pleurisy (23.62%). Concurrent PTB-EPTB was found in 129,422 cases (29.48%) of which tuberculous pleurisy was the most common concurrent EPTB type observed. The multivariable logistic regression models demonstrated that odds ratios of concurrent PTB-EPTB cases varied by gender and age group. For PTB cases with concurrent EPTB, the strongest association was found between PTB and concurrent bronchial tuberculosis (lift = 1.09). For EPTB cases with concurrent PTB, the strongest association was found between pharyngeal/laryngeal tuberculosis and concurrent PTB (lift = 1.11). Confidence and lift values of concurrent PTB-EPTB cases varied with gender and age. Conclusions: Numerous concurrent PTB-EPTB case types were observed, with confidence and lift values varying with gender and age. Clinicians should screen for concurrent PTB-EPTB in order to improve treatment outcomes.

Epidemiology of concurrent extrapulmonary tuberculosis in inpatients with extrapulmonary tuberculosis lesions in China: a large-scale observational multi-centre investigation.

AIMS: A high proportion of all patients with tuberculosis (TB) present with extrapulmonary TB (EPTB), including concurrent EPTB involving more than one extrapulmonary lesion site. However, previous reports only characterized lesions of single-site EPTB cases. This study aimed to investigate epidemiological characteristics and association rules of concurrent EPTB cases in China. METHODS: An observational multi-centre study of 208,214 patients with EPTB lesions was undertaken in China from January 2011 to December 2017. Multi-variable logistic regression analysis was used to identify associations between gender and concurrent EPTB, and age and concurrent EPTB. Association rules were analysed for significance using the Apriori algorithm. RESULTS: The most common EPTB lesion was tuberculous pleurisy (49.8%), followed by bronchial TB (14.8%) and tuberculous meningitis (7.6%). The most common type of concurrent EPTB was tuberculous pleurisy concurrent with tuberculous peritonitis (1.80%). In total, 22 association rules, including 20 strong association rules, were identified; among these, the highest confidence rates were found for tuberculous myelitis concurrent with tuberculous meningitis, and sacral TB concurrent with lumbar vertebral TB. The association rules of EPTB concurrent with other EPTB types were found to vary with gender and age. The confidence rate of tuberculous myelitis concurrent with tuberculous meningitis was higher in females (83.67%) than males, and was highest in patients aged 25-34 years (87.50%). CONCLUSIONS: Many types of concurrent EPTB were found. Greater awareness of concurrent EPTB disease characteristics is needed to ensure timely clinical diagnosis and treatment of this disease.

The prevalence and risks of major comorbidities among inpatients with pulmonary tuberculosis in China from a gender and age perspective: a large-scale multicenter observational study.

In clinical practice, PTB patients have concurrent many types of comorbidities such as pneumonia, liver disorder, diabetes mellitus, hematological disorder, and malnutrition. Detecting and treating specific comorbidities and preventing their development are important for PTB patients. However, the prevalence of most comorbid conditions in patients with PTB is not well described. We conducted a large-scale, multicenter, observational study to elucidate and illustrate the prevalence rates of major comorbidities in inpatients at 21 hospitals in China. The 19 specific comorbidities were selected for analysis in this patient cohort, and stratified the inpatient cohort according to age and gender. A total of 355,929 PTB inpatients were included, with a male:female ratio of 1.98 and the proportion of ≥ 65 years PTB inpatients was the most. Approximately 70% of PTB inpatients had at least one defined type of comorbidity. The prevalence of 19 specific comorbidities in inpatients with PTB was analyzed, with pneumonia being the most common comorbidity. The prevalence of most comorbidities was higher in males with PTB except thyroid disorders, mental health disorders, etc. The prevalence of defined most comorbidities in patients with PTB tended to increase with increasing age, although some specific comorbidities tended to increase initially then decrease with increasing age. Our study describes multiple clinically important comorbidities among PTB inpatients, and their prevalence between different gender and age groups. The results will enhance the clinical aptitude of physicians who treat patients with PTB to recognize, diagnose, and treat PTB comorbidities early.

The epidemiology of extrapulmonary tuberculosis in China: A large-scale multi-center observational study.

Tuberculosis (TB) remains a serious global public health problem in the present. TB also affects other sites (extrapulmonary tuberculosis, EPTB), and accounts for a significant proportion of tuberculosis cases worldwide. In order to comprehensively understand epidemiology of EBTB in China, and improve early diagnosis and treatment, we conducted a large-scale multi-center observational study to assess the demographic data and the prevalence of common EPTB inpatients, and further evaluate the prevalence of EPTB concurrent with Pulmonary tuberculosis (PTB) and the associations between multiple EPTB types and gender-age group in China. All consecutive age≥15yr inpatients with a confirmed diagnosis of EPTB during the period from January 2011 to December 2017 were included in the study. The descriptive statistical analysis included median and quartile measurements for continuous variables, and frequencies and proportions with 95% confidence intervals (CIs) for categorical variables. Multinomial logistic regression analysis was used to compare the association of multiple EPTB types between age group and gender. The results showed that the proportion of 15-24 years and 25-34 years in EPTB inpatients were the most and the ratio of male: female was 1.51. Approximately 70% of EPTB inpatients were concurrent with PTB or other types of EPTB. The most common of EPTB was tuberculous pleurisy (50.15%), followed by bronchial tuberculosis (14.96%), tuberculous lymphadenitis of the neck (7.24%), tuberculous meningitis (7.23%), etc. It was found that many EPTB inpatients concurrent with PTB. The highest prevalence of EPTB concurrent with PTB was pharyngeal/laryngeal tuberculosis (91.31%), followed by bronchial tuberculosis (89.52%), tuberculosis of hilar lymph nodes (79.52%), tuberculosis of mediastinal lymph nodes (79.13%), intestinal tuberculosis (72.04%), tuberculous pleurisy (65.31%) and tuberculous meningitis (62.64%), etc. The results from EPTB concurrent with PTB suggested that females EPTB inpatients were less likely to be at higher risk of concurrent PTB (aOR = 0.819, 95%CI:0.803-0.835) after adjusted by age. As age increasing, the trend risk of concurrent PTB decreased (aOR = 0.994, 95%CI: 0.989-0.999) after adjusted by gender. Our study demonstrated that the common EPTB were tuberculous pleurisy, bronchial tuberculosis, tuberculous lymphadenitis of the neck, tuberculous meningitis, etc. A majority of patients with pharyngeal/laryngeal tuberculosis, bronchial tuberculosis, tuberculosis of hilar/mediastinal lymph nodes, intestinal tuberculosis, tuberculous pleurisy, tuberculous meningitis, etc. were concurrent with PTB. Female EPTB inpatients were less likely to be at higher risk of concurrent PTB, and as age increasing, the trend risk of concurrent PTB decreased. The clinicians should be alert to the presence of concurrent tuberculosis in EPTB, and all suspected cases of EPTB should be assessed for concomitant PTB to determine whether the case is infectious and to help for early diagnosis and treatment.

Mössbauer properties of the diferric cluster and the differential iron(II)-binding affinity of the iron sites in protein R2 of class Ia Escherichia coli ribonucleotide reductase: a DFT/electrostatics study.

The R2 subunit of class-Ia ribonucleotide reductase (RNR) from Escherichia coli (E. coli) contains a diiron active site. Starting from the apo-protein and Fe(II) in solution at low Fe(II)/apoR2 ratios, mononuclear Fe(II) binding is observed indicating possible different Fe(II) binding affinities for the two alternative sites. Further, based on their Mössbauer spectroscopy and two-iron-isotope reaction experiments, Bollinger et al. (J. Am. Chem. Soc., 1997, 119, 5976-5977) proposed that the site Fe1, which bonds to Asp84, should be associated with the higher observed (57)Fe Mössbauer quadrupole splitting (2.41 mm s(-1)) and lower isomer shift (0.45 mm s(-1)) in the Fe(III)Fe(III) state, site Fe2, which is further from Tyr122, should have a greater affinity for Fe(II) binding than site Fe1, and Fe(IV) in the intermediate X state should reside at site Fe2. In this paper, using density functional theory (DFT) incorporated with the conductor-like screening (COSMO) solvation model and with the finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) methodologies, we have demonstrated that the observed large quadrupole splitting for the diferric state R2 does come from site Fe1(III) and it is mainly caused by the binding position of the carboxylate group of the Asp84 sidechain. Further, a series of active site clusters with mononuclear Fe(II) binding at either site Fe1 or Fe2 have been studied, which show that with a single dielectric medium outside the active site quantum region, there is no energetic preference for Fe(II) binding at one site over another. However, when including the explicit extended protein environment in the PB-SCRF model, the reaction field favors the Fe(II) binding at site Fe2 rather than at site Fe1 by ~9 kcal mol(-1). Therefore our calculations support the proposal of the previous Mössbauer spectroscopy and two-iron-isotope reaction experiments by Bollinger et al.

Geometric and electrostatic study of the [4Fe-4S] cluster of adenosine-5'-phosphosulfate reductase from broken symmetry density functional calculations and extended X-ray absorption fine structure spectroscopy.

Adenosine-5'-phosphosulfate reductase (APSR) is an iron-sulfur protein that catalyzes the reduction of adenosine-5'-phosphosulfate (APS) to sulfite. APSR coordinates to a [4Fe-4S] cluster via a conserved CC-X(~80)-CXXC motif, and the cluster is essential for catalysis. Despite extensive functional, structural, and spectroscopic studies, the exact role of the iron-sulfur cluster in APS reduction remains unknown. To gain an understanding into the role of the cluster, density functional theory (DFT) analysis and extended X-ray fine structure spectroscopy (EXAFS) have been performed to reveal insights into the coordination, geometry, and electrostatics of the [4Fe-4S] cluster. X-ray absorption near-edge structure (XANES) data confirms that the cluster is in the [4Fe-4S](2+) state in both native and substrate-bound APSR while EXAFS data recorded at ~0.1 Å resolution indicates that there is no significant change in the structure of the [4Fe-4S] cluster between the native and substrate-bound forms of the protein. On the other hand, DFT calculations provide an insight into the subtle differences between the geometry of the cluster in the native and APS-bound forms of APSR. A comparison between models with and without the tandem cysteine pair coordination of the cluster suggests a role for the unique coordination in facilitating a compact geometric structure and "fine-tuning" the electronic structure to prevent reduction of the cluster. Further, calculations using models in which residue Lys144 is mutated to Ala confirm the finding that Lys144 serves as a crucial link in the interactions involving the [4Fe-4S] cluster and APS.

DFT calculations for intermediate and active states of the diiron center with a tryptophan or tyrosine radical in Escherichia coli ribonucleotide reductase.

Class Ia ribonucleotide reductase subunit R2 contains a diiron active site. In this paper, active-site models for the intermediate X-Trp48(•+) and X-Tyr122(•), the active Fe(III)Fe(III)-Tyr122(•), and the met Fe(III)Fe(III) states of Escherichia coli R2 are studied, using broken-symmetry density functional theory incorporated with the conductor-like screening solvation model. Different structural isomers and different protonation states have been explored. Calculated geometric, energetic, Mössbauer, hyperfine, and redox properties are compared with available experimental data. Feasible detailed structures of these intermediate and active states are proposed. Asp84 and Trp48 are most likely the main contributing residues to the result that the transient Fe(IV)Fe(IV) state is not observed in wild-type class Ia E. coli R2. Asp84 is proposed to serve as a proton-transfer conduit between the diiron cluster and Tyr122 in both the tyrosine radical activation pathway and the first steps of the catalytic proton-coupled electron-transfer pathway. Proton-coupled and simple redox potential calculations show that the kinetic control of proton transfer to Tyr122(•) plays a critical role in preventing reduction from the active Fe(III)Fe(III)-Tyr122(•) state to the met state, which is potentially the reason why Tyr122(•) in the active state can be stable over a very long period.

Density functional theory analysis of structure, energetics, and spectroscopy for the Mn-Fe active site of Chlamydia trachomatis ribonucleotide reductase in four oxidation states.

Models for the Mn-Fe active site structure of ribonucleotide reductase (RNR) from pathogenic bacteria Chlamydia trachomatis (Ct) in different oxidation states have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and also with finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. The detailed structures for the reduced Mn(II)-Fe(II), the met Mn(III)-Fe(III), the oxidized Mn(IV)-Fe(III) and the superoxidized Mn(IV)-Fe(IV) states are predicted. The calculated properties, including geometries, (57)Fe Mossbauer isomer shifts and quadrupole splittings, and (57)Fe and (55)Mn electron nuclear double resonance (ENDOR) hyperfine coupling constants, are compared with the available experimental data. The Mössbauer and energetic calculations show that the (mu-oxo, mu-hydroxo) models better represent the structure of the Mn(IV)-Fe(III) state than the di-mu-oxo models. The predicted Mn(IV)-Fe(III) distances (2.95 and 2.98 A) in the (mu-oxo, mu-hydroxo) models are in agreement with the extended X-ray absorption fine structure (EXAFS) experimental value of 2.92 A (Younker et al. J. Am. Chem. Soc. 2008, 130, 15022-15027). The effect of the protein and solvent environment on the assignment of the Mn metal position is examined by comparing the relative energies of alternative mono-Mn(II) active site structures. It is proposed that if the Mn(II)-Fe(II) protein is prepared with prior addition of Mn(II) or with Mn(II) richer than Fe(II), Mn is likely positioned at metal site 2, which is further from Phe127.

Quantum cluster size and solvent polarity effects on the geometries and Mössbauer properties of the active site model for ribonucleotide reductase intermediate X: a density functional theory study.

In studying the properties of metalloproteins using ab initio quantum mechanical methods, one has to focus on the calculations on the active site. The bulk protein and solvent environment is often neglected, or is treated as a continuum dielectric medium with a certain dielectric constant. The size of the quantum cluster of the active site chosen for calculations can vary by including only the first-shell ligands which are directly bound to the metal centers, or including also the second-shell residues which are adjacent to and normally have H-bonding interactions with the first-shell ligands, or by including also further hydrogen bonding residues. It is not well understood how the size of the quantum cluster and the value of the dielectric constant chosen for the calculations will influence the calculated properties. In this paper, we have studied three models (A, B, and C) of different sizes for the active site of the ribonucleotide reductase intermediate X, using density functional theory (DFT) OPBE functional with broken-symmetry methodology. Each model is studied in gas-phase and in the conductor-like screening (COSMO) solvation model with different dielectric constants ε = 4, 10, 20, and 80, respectively. All the calculated Fe-ligand geometries, Heisenberg J coupling constants, and the Mössbauer isomer shifts, quadrupole splittings, and the (57)Fe, (1)H, and (17)O hyperfine tensors are compared. We find that the calculated isomer shifts are very stable. They are virtually unchanged with respect to the size of the cluster and the dielectric constant of the environment. On the other hand, certain Fe-ligand distances are sensitive to both the size of the cluster and the value of ε. ε = 4, which is normally used for the protein environment, appears too small when studying the diiron active site geometry with only the first-shell ligands as seen by comparisons with larger models.

DFT calculations of comparative energetics and ENDOR/Mössbauer properties for two protonation states of the iron dimer cluster of ribonucleotide reductase intermediate X.

Two models (I and II) for the active site structure of class-I ribonucleotide reductase (RNR) intermediate X in subunit R2 have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and with the finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. Only one of the bridging groups between the two iron centers is different between model-I and model-II. Model-I contains two mu-oxo bridges, while model-II has one bridging oxo and one bridging hydroxo. These are large active site models including up to the fourth coordination shell H-bonding residues. Mössbauer and ENDOR hyperfine property calculations show that model-I is more likely to represent the active site structure of RNR-X. However, energetically our pK(a) calculations at first highly favored the bridging oxo and hydroxo (in model-II) structure of the diiron center rather than having the di-oxo bridge (in model-I). Since the Arg236 and the nearby Lys42, which are very close to the diiron center, are on the protein surface of RNR-R2, it is highly feasible that one or two anion groups in solution would interact with the positively charged side chains of Arg236 and Lys42. The anion group(s) can be a reductant, phosphate, sulfate, nitrate, and other negatively charged groups existing in biological environments or in the buffer of the experiment. Since sulfate ions certainly exist in the buffer of the ENDOR experiment, we have examined the effect of the sulfate (SO(4)(2-), surrounded by explicit water molecules) H-bonding to the side chain of Arg236. We find that when sulfate interacts with Arg236, the carboxylate group of Asp237 tends to be protonated, and once Asp237 is protonated, the Fe(iii)Fe(iv) center in X favors the di-oxo bridge (model-I). This would explain that the ENDOR observed RNR-X active site structure is likely to be represented by model-I rather than model-II.

Structural model studies for the peroxo intermediate P and the reaction pathway from P-->Q of methane monooxygenase using broken-symmetry density functional calculations.

Several structural models for the active site of the peroxo intermediate state "P" of the hydroxylase component of soluble methane monooxygenase (MMOH) have been studied, using two DFT functionals OPBE and PW91 with broken-symmetry methodology and the conductor-like screening (COSMO) solvation model. These active site models have different O2 binding modes to the diiron center, such as the mu-eta2,eta2, trans-mu-1,2 and cis-mu-1,2 conformations. The calculated properties, including optimized geometries, electronic energies, Fe net spin populations, and Mössbauer isomer shift and quadrupole splitting values, have been reported and compared with available experimental results. The high-spin antiferromagnetically (AF) coupled Fe3+ sites are correctly predicted by both OPBE and PW91 methods for all active site models. Our data analysis and comparisons favor a cis-mu-1,2 structure (model cis-mu-1,2a shown in Figure 9) likely to represent the active site of MMOH-P. Feasible structural changes from MMOH-P to another intermediate state MMOH-Q are also proposed, where the carboxylate group of Glu243 side chain has to open up from the mono-oxygen bridging position, and the dissociations of the terminal H2O ligand from Fe1 and of the oxygen atom in the carboxylate group of Glu144 from Fe2 are also necessary for the O2 binding mode changes from cis to trans. The O-O bond is proposed to break in the trans-conformation and forms two mu-oxo bridges in MMOH-Q. The terminal H2O molecule and the Glu144 side chain then rebind with Fe1 and Fe2, respectively, in Q.

Structural Model Studies for the High-Valent Intermediate Q of Methane Monooxygenase from Broken-Symmetry Density Functional Calculations.

Mössbauer isomer shift parameters have been obtained for both density functional theory (DFT) OPBE and OLYP functionals by linear regressions between the measured isomer shifts and calculated electron densities at Fe nuclei for a number of Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes grouped separately. The calculated isomer shifts and quadrupole splittings on the sample Fe complexes from OPBE and OLYP functionals are similar to those of PW91 calculations (J. Comput. Chem. 27 (2006) 1292), however the fit parameters from the linear regressions differ between PW91 and OPBE, OLYP. Four models for the active site structure of the hydroxylase component of soluble methane monooxygenase (MMOH) have been studied, using three DFT functionals OPBE, OLYP, and PW91, incorporated with broken-symmetry methodology and the conductor-like screening (COSMO) solvation model. The calculated properties, including optimized geometries, electronic energies, pK(a)'s, Fe net spin populations, and Mössbauer isomer shifts and quadrupole splittings, have been reported and compared with available experimental values. The high-spin antiferromagnetically (AF) coupled Fe(4+) sites are correctly predicted by OPBE and OLYP methods for all active site models. PW91 potential overestimates the Fe-ligand covalencies for some of the models because of spin crossover. Our calculations and data analysis support the structure (our current model II shown in Figure 8) proposed by Friesner and Lippard's group (J. Am. Chem. Soc. 123 (2001) 3836-3837), which contains an Fe(4+)(µ-O)(2)Fe(4+) center, one axial water which also H-bonds to both side chains of Glu243 and Glu114, and one bidentate carboxylate group from the side chain of Glu144, which is likely to represent the active site of MMOH-Q. A new model structure (model IV shown in Figure 9), which has a terminal hydroxo and a protonated His147 which is dissociated from a nearby Fe, is more asymmetric in its Fe(µ-O)(2)Fe diamond core, and is another very good candidate for intermediate Q.

Experimental and DFT studies: novel structural modifications greatly enhance the solvent sensitivity of live cell imaging dyes.

Structural modifications of previously reported merocyanine dyes (Toutchkine, A.; Kraynov, V.; Hahn, K. J. Am. Chem. Soc. 2003, 125, 4132-4145) were found to greatly enhance the solvent dependence of their absorbance and fluorescence emission maxima. Density functional theory (DFT) calculations have been performed to understand the differences in optical properties between the new and previously synthesized dyes. Absorption and emission energies were calculated for several new dyes using DFT vertical self-consistent reaction field (VSCRF) methods. Geometries of ground and excited states were optimized with a conductor-like screening model (COSMO) and self-consistent field (SCF) methods. The new dyes have enhanced zwitterionic character in the ground state and much lower polarity in the excited state, as shown by the DFT-VSCRF calculations. Consistently, the position of the absorption bands are strongly blue-shifted in more polar solvent (methanol compared to benzene), as predicted by the DFT spectral calculations. Inclusion of explicit H-bonding solvent molecules within the quantum model further enhances the predicted shifts and is consistent with the observed spectral broadening. Smaller but significant spectral shifts in polar versus nonpolar solvent are predicted and observed for emission bands. The new dyes show large fluorescence quantum yields in polar hydrogen-bonding solvents; qualitatively, the longest bonds along the conjugated chain at the excited S1 state minimum are shorter in the more polar solvent, inhibiting photoisomerization. The loss of photostability of the dyes is a consequence of the reaction with and electron transfer to singlet oxygen, starting oxidative dye cleavage. The calculated vertical ionization potentials of three dyes I-SO, AI-SO(4), and AI-BA(4) in benzene and methanol are consistent with their relative photobleaching rates; the charge distributions along the conjugated chains for the three dyes are similarly predictive of higher reaction rates for AI-SO(4) and AI-BA(4) than for I-SO. Time-dependent DFT calculations were also performed on AI-BA(4); these were less accurate than the VSCRF method in predicting the absorption energy shift from benzene to methanol.

Density functional theory study of Fe(IV) d-d optical transitions in active-site models of class I ribonucleotide reductase intermediate X with vertical self-consistent reaction field methods.

The Fe(IV) d-d transition energies for four active-site structural models of class I ribonucleotide reductase (RNR) intermediate X have been calculated using broken-symmetry density functional theory incorporated with the Slater transition state vertical self-consistent reaction field methodology. Our model I (Figure 1), which contains two mu-oxo bridges, one terminal water, and one bidentate carboxylate group, yields the best Fe(IV) d-d transition energies compared with experiment. Our previous study (J. Am. Chem. Soc. 2005, 127, 15778-15790) also shows that most of the other calculated properties of model I in both native and mutant Y122F forms, including geometries, spin states, pKa's, 57Fe, 1H, and 17O hyperfine tensors, and 57Fe Mössbauer isomer shifts and quadrupole splittings, are also the best in agreement with the available experimental data. This model is likely to represent the active-site structure of the intermediate state X of RNR.

Structure, redox, pKa, spin. A golden tetrad for understanding metalloenzyme energetics and reaction pathways.

After a review of the current status of density functional theory (DFT) for spin-polarized and spin-coupled systems, we focus on the resting states and intermediates of redox-active metalloenzymes and electron transfer proteins, showing how comparisons of DFT-calculated spectroscopic parameters with experiment and evaluation of related energies and geometries provide important information. The topics we examine include (1) models for the active-site structure of methane monooxygenase intermediate Q and ribonucleotide reductase intermediate X; (2) the coupling of electron transfer to proton transfer in manganese superoxide dismutase, with implications for reaction kinetics; (3) redox, pK(a), and electronic structure issues in the Rieske iron-sulfur protein, including their connection to coupled electron/proton transfer, and an analysis of how partial electron delocalization strongly alters the electron paramagnetic resonance spectrum; (4) the connection between protein-induced structural distortion and the electronic structure of oxidized high-potential 4Fe4S proteins with implications for cluster reactivity; (5) an analysis of cluster assembly and central-atom insertion into the FeMo cofactor center of nitrogenase based on DFT structural and redox potential calculations.

DFT calculations of 57Fe Mössbauer isomer shifts and quadrupole splittings for iron complexes in polar dielectric media: applications to methane monooxygenase and ribonucleotide reductase.

To predict the isomer shifts of Fe complexes in different oxidation and spin states more accurately, we have performed linear regression between the measured isomer shifts (delta(exp)) and DFT (PW91 potential with all-electron triple-zeta plus polarization basis sets) calculated electron densities at Fe nuclei [rho(0)] for the Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes separately. The geometries and electronic structures of all complexes in the training sets are optimized within the conductor like screening (COSMO) solvation model. Based on the linear correlation equation delta(exp) = alpha[rho(0) - 11884.0] + C, the best fitting for 17 Fe(2+,2.5+) complexes (totally 31 Fe sites) yields alpha = -0.405 +/- 0.042 and C = 0.735 +/- 0.047 mm s(-1). The correlation coefficient is r = -0.876 with a standard deviation of SD = 0.075 mm s(-1). In contrast, the linear fitting for 19 Fe(2.5+,3+,3.5+,4+) complexes (totally 30 Fe sites) yields alpha = -0.393 +/- 0.030 and C = 0.435 +/- 0.014 mm s(-1), with the correlation coefficient r = -0.929 and a standard deviation SD = 0.077 mm s(-1). We provide a physical rationale for separating the Fe(2+,2.5+) fit from the Fe(2.5+,3+,3.5+,4+) fit, which also is clearly justified on a statistical empirical basis. Quadrupole splittings have also been calculated for these systems. The correlation between the calculated (DeltaE(Q(cal))) and experimental (DeltaE(Q(exp))) quadrupole splittings based on |DeltaE(Q(exp))| = A |DeltaE(Q(cal))| + B yields slope A, which is almost the ideal value 1.0 (A = 1.002 +/- 0.030) and intercept B almost zero (B = 0.033 +/- 0.068 mm s(-1)). Further calculations on the reduced diferrous and oxidized diferric active sites of class-I ribonucleotide reductase (RNR) and the hydroxylase component of methane monooxygenase (MMOH), and on a mixed-valent [(tpb)Fe3+(mu-O)(mu-CH3CO2)Fe4+(Me3[9]aneN3)]2+ (S = 3/2) complex and its corresponding diferric state have been performed. Calculated results are in very good agreement with the experimental data.

Seven clues to the origin and structure of class-I ribonucleotide reductase intermediate X.

Class-I ribonucleotide reductases (RNRs) are aerobic enzymes that catalyze the reduction of ribonucleotides to deoxyribonucleotides providing the required building blocks for DNA replication and repair. These ribonucleotide-to-deoxyribonucleotide reactions occur by a long range radical (or proton-coupled-electron-transfer) propagation mechanism initiated by a fairly stable tyrosine radical ("the pilot light"). When this pilot light goes out, the tyrosine radical is regenerated by a high-oxidation-state enzyme intermediate, called X. The active site of class-I RNR-X has been recognized as a spin coupled Fe(III)Fe(IV) center with S(total)=1/2 ground state. Although several clues have been obtained from Mössbauer, (57)Fe, (1)H, (17)O(2), and H(2)(17)O ENDOR (electron-nuclear double resonance), EXAFS (extended X-ray absorption fine structure), and MCD (magnetic circular dichroism) experiments, the detailed structure of the intermediate X is still unknown. In the past three years, we have been studying the properties of a set of model clusters for RNR-X using broken-symmetry density functional theory (DFT), and have compared them with the available experimental results. Based on the detailed analysis and comparisons, we have proposed a definite form for the active site structure of class-I RNR intermediate X. The puzzle is now set: can you find any flaws in the argument or evidence? Can you add anything further to the current experimental picture? The argument is formulated from seven experimental clues with associated calculations and models.

Active site structure of class I ribonucleotide reductase intermediate X: a density functional theory analysis of structure, energetics, and spectroscopy.

Several models for the active site structure of class I ribonucleotide reductase (RNR) intermediate X have been studied in the work described in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor-like screening (COSMO) solvation model. The calculated properties, including geometries, spin states, 57Fe, 1H, and 17O hyperfine tensors, Mössbauer isomer shifts, and quadrupole splittings, and the estimation of the Fe(IV) d-d transition energies have been compared with the available experimental values. On the basis of the detailed analysis and comparisons, we propose a definite form for the active site structure of class I RNR intermediate X, which contains an Fe1(III)Fe2(IV) center (where Fe1 is the iron site closer to Tyr122, and the two iron sites are high-spin antiferromagnetically coupled to give a total 1/2 net spin), two mu-oxo bridges, one terminal water which binds to Fe1(III) and also H-bonds to both side chains of Asp84 and Glu238, and one bidentate carboxylate group from the side chain of Glu115.